Effect of Inhomogeneity on s-wave Superconductivity in the Attractive Hubbard 

Model 
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Inhomogeneous s-wave superconductivity is studied in the two-dimensional, square lattice attrac- 
tive Hubbard Hamiltonian using the Bogoliubov-de Gennes (BdG) mean field approximation. We 
find that at weak coupling, and for densities mainly below half-filling, an inhomogeneous interaction 
in which the on-site interaction Ui takes on two values, Ui = 0, 2(7 results in a larger zero tempera- 
ture pairing amplitude, and that the superconducting Tc can also be significantly increased, relative 
to a uniform system with Ui = U on all sites. These effects are observed for stripe, checkerboard, 
and even random patterns of the attractive centers, suggesting that the pattern of inhomogeneity is 
unimportant. Monte Carlo calculations which reintroduce some of the fluctuations neglected within 
the BdG approach see the same effect, both for the attractive Hubbard model and a Hamiltonian 
with d-wave pairing symmetry. 

PACS numbers: 



I. MOTIVATION 



One of the main themes of recent studies of strongly 
correlated electronic systems is the importance of spatial 
inhomogeneities. These can result either from intrinsic 
quenched disorder in the system, as in the metal-insulator 
transition in two dimensionsli^, or arise spontaneously in 
an otherwise translationally invariant system. For ex- 
ample, holes doped into the high temperature supercon- 
ductors, (HTS) appear not to spread out uniformly in 
the Cu02 planes, but instead arrange themselves in the 
form of stripes, checkerboard or perhaps even more ex- 
otic structures i^i^i^iSii Besides the cuprate superconduc- 
tors, such spatially varying density and spin structures 
are also key features in the physics of the manganiteaSiS 
and cobaltitesii2iii 

Considerable theoretical work on the interplay be- 
tween spatial inhomogeneity, magnetism, and supercon- 
ductivity has utilized the repulsive Hubbard and t-J 
Hamiltoniansi2iiii^iiSiiSiiLia For the 2D square lattice 
these models are known to display antiferromagnetism at 
half-filling, and, although it is less certain, perhaps also 
d-wave superconductivity when doped. There is consid- 
erable evidence that they also might possess inhomoge- 
neous stripe or checkerboard ground states>i^*i& While 
DMRG treatmentsi^ provide detailed information on the 
real space charge, spin, and pairing orders, the precise 
nature of the interplay, and whether the different orders 
compete or cooperate, remains unclear. In addition, the 
enhancement of the superconducting transition temper- 
ature Tc by local inhomogeneity has been demonstrated 
by Martin et ali^ 



II. MODEL AND METHODOLOGY 

In this paper we address the general issue of whether 
inhomogeneous regions of attraction favor superconduc- 
tivity relative to the homogeneous system with the same 
average attraction, either by increasing the zero temper- 
ature pairing amplitude or the transition temperature. 
In many of the systems for which this question is funda- 
mental, such as the cuprate superconductors mentioned 
above, the situation is complicated by the presence of 
other types of order such as antiferromagnetism, exotic 
spin-gap phases, and nontrivial d-wave symmetry of the 
superconducting order parameter. Rather than using a 
model like the repulsive Hubbard Hamiltonian which in- 
corporates this full richness, it is useful to study the prob- 
lem first in a more simple context. Here we will present 
a solution of the inhomogeneous Bogoliubov-de Gennes 
(BdG) equations for the attractive Hubbard Hamilto- 
nian, 

<ij>,(7 

- ^^^4aC■^,<y -^\U^\n^]n,^ , (1) 

ia i 

with t the hopping amplitude, /i the chemical poten- 
tial and Ui the local attractive interaction between the 
fermions of opposite spins residing on the same lattice site 
i. Our focus will be on inhomogeneous patterns in the in- 
teraction Ui. The interaction in the attractive Hubbard 
model can be thought of as a phenomenological one, origi- 
nating, for example, from integrating out a local phonon 
modc.^*^ The two-dimensional uniform attractive Hub- 
bard model is known to yield degenerate superconductiv- 
ity and charge density wave (CDW) long range order at 
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half-filling and zero temperature i2ii2242^ However, away 
from half-filling, the CDW pairing symmetry is broken 
and superconductivity is more favorable, and the super- 
conducting phase transition is at finite temperature. 

Within the BdG mean field decomposition, we replace 
the local pairing amplitude and local density by their 
average values, — (^CiiCii) and (njo-) — (cJo-Cio-) and 
arrive at the quadratic effective Hamiltonian 



n 



< ij > , cr ia 

^'U,\[A4^cl + A*c,^c,^], (2) 



Ui\{ni)/2 includes a site-dependent 
[^i) = Y.ai'^i-y)- energies wiU be 



where fii = ji + 
Hartree shift with 
referenced to < = 1. 

We adopt the criterion of comparing the tendency for 
superconductivity in the homogeneous system with the 
same attraction — [/ on all lattice sites, with cases when 
sites with attraction are mixed with sites where the at- 
traction is absent, i.e., Ui = 0ii2i24 Specifically, we have 
studied systems in which sites with attractive interaction 
are randomly distributed^ or arranged in checkerboard 
and stripe patterns. In all three inhomogencous patterns, 
exactly half of the lattice sites carry interaction, and the 
interacting sites carry twice the value of U as in the case 
of the uniform pattern. 

This conventional mean-field approach does not cap- 
ture the Kosterlitz-Thouless nature of the phase transi- 
tion in two dimensions. Nevertheless, this weakness can 
be repaired^ upon regarding the local pairing amplitudes 
as complex variables and performing a finite temperature 
Monte Carlo integration over the associated amplitude 
and phase degrees of freedom. Unlike BCS, this Monte 
Carlo mean field (MCMF) approach, allows identification 
of the weak and strong coupling regimes via the phase 
correlation function. We will use this Monte Carlo tech- 
nique as an independent confirmation of our results. 



III. SOLUTIONS AND PHASE DIAGRAM 

For the doping of n = 0.8 at T = 0, as depicted 
in Fig. n^a), an inhomogencous system with bimodal 
Ui = 0, 2U has a larger zero temperature gap than a uni- 
form system with Ui = U below the maximum value of 
t/(r=i) ~ 2.0 with r = Ai„homog/Auniform- At half-fiUiug 
in Fig. njb),this effect for the stripe pattern terminates 
at a slightly smaller value of C^(r=i) ~ 1-5 and at the same 
time, the magnitude of A^j.^]^) has also diminished. There 
is an upper critical Uc above which superconductivity is 
obliterated by inhomogeneity. For the checkerboard pat- 
tern at half-filling in particular, A becomes zero as one 
enters the insulating CDW phase of static pairs. For the 
doping of n = 1.2 (above half-filling) in Fig. ^c) for 
random and stripe patterns and A for the checkerboard 
pattern still remains infinitesimal over almost the entire 
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FIG. 1: Variation of the mean zero temperature pairing am- 
plitude A, versus mean interaction [7, for (a) below half-filling 
(n = 0.8), (b) half-filling (n = 1.0), and (c) above half-filling 
(n = 1.2). There is no pairing at n > 1 in the checker- 
board case. For the stripe, checkerboard, and uniform pat- 
terns, data are fully converged on the momentum space grid. 
That is, data for the indicated lattice sizes lie on top of each 
other. For the random pattern it is not possible to study 
large lattices, and a single size, L = 14 is shown. Panel (d) 
compares the MCMF results (T 0) at L = 10 with their 
BdG counterparts at L = 10, 500 and 1500 for the uniform 
and checkerboard patterns. 



range of U . The MCMF data in Fig. ^d), provide con- 
firmation of the BdG results. While the MCMF method 
remains limited to relatively small lattice sizes, BdG al- 
lows us to invoke lattice sizes as large as L = 1500 to re- 
duce the finite size effects at small U values. As observed 
in Fig. n^d) , the finite size effects are more significant for 
the uniform pattern when [/ < 3.0 and less severe for the 
checkerboard. The good agreement between these two 
approaches helps justify both the application and results 
of the BdG technique. 

Fig. 121 illustrates further the size of the effect due to 
inhomogeneity, using data over a broad range of densi- 
ties. For the checkerboard pattern in Fig. [2Ia), as the 
doping changes from n = 0.2 to n = 0.9, the intersection 
point of the ratio r — Achcck/Auniform and unity line is 
shifted towards smaller U(^^i)- Precisely at half-filling, 
A for the checkerboard pattern vanishes. For stripes, 
however, as shown in Fig. |2Ib), the enhancement shift 
towards smaller C/(i.=i]_ continues through half-filling to 
almost n = 1.35. Fig. |2jc) shows the value of the pairing 
amplitude A(j.^x') as a function of the electron doping. 
Fillings around quarter-filling (n = 0.5) have the largest 
ratio in pair amplitude of the bimodal Ui = 0, 2U to 
uniform Ui — U interaction distributions. 

The lack of superconductivity in the checkerboard pat- 
tern at half-filling is associated with the formation of a 
competing, insulating CDW phase. This can be seen di- 
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FIG. 2: Panel (a): the ratio r — Achcck/Auniform as a 
function of U and different dopings of electrons at T = 0. 
Panel (b): the same results as for panel (a) for the ratio 
r = Astripc/Aunitorm- Panel (c): Value of the pairing am- 
plitude at which r — 1 versus the electron doping. At weak 
coupling, the zero temperature pairing amplitude can be a 
factor of three to four larger than the uniform system. Lat- 
tice sizes of L = 1500 for the uniform and L = 500 and 1000 
for inhomogeneous patterns were utilized. 
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FIG. 3: Variation of the pairing amplitude versus the mean 
interaction U at T = 0. Panel (a) shows A as the average 
of Au and A(u=o) for the checkerboard pattern at n = 0.8 
and the same for stripes in panel (b). The enhancement of 
pairing by inhomogeneity is seen to arise from a proximity 
effect whereby significant pair amplitude is induced on the 
U = sites by the U ^ sites. Lattice sizes of L = 500 and 
1000 were utilized. 
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FIG. 4: Variation of the pairing amplitude versus temperature 
T for n = 0.8. Panel (a): A as a function of T for [7 = 1.5 
where checkerboard and stripe patterns with Ui = 0, 2U have 
a larger T — pair amplitude than Ui = (Fig. Q^a)). Panel 
(b): the same results for U = 3.0. Lattice sizes of L = 1500 
for the uniform and L — 500 and 1000 for inhomogeneous 
patterns were utilized. 



rectly in the density of states as a gap develops at the 
Fermi energy despite the fact that Age = 0. The occupa- 
tion of sites in real space becomes increasingly disparate 
as U increases, with the sites with U becoming fuUy 
packed with ny « 2.0, while the non-interacting sites be- 
come empty, n(u=o) ~ 0.0. For stripes at half-filhng also, 
at large enough values of U superconductivity is obliter- 
ated. Above half-filling, however, when the order pa- 
rameter vanishes for both the checkerboard and stripes, 
the density of states remains finite at the Fermi energy, 
indicating a metallic phase. 

The proximity effect for the non-interacting sites 
neighbored by the interacting sites plays a major role 
in the magnitude of the pair amplitude in the inho- 
mogeneous lattice. In Fig. ^a), {n = 0.8) the value 
of U at which the checkerboard and stripe patterns 
have the same pair amplitude as the uniform system is 
t^(r=i) ~ 2.0. In Fig. ISl^a) and (b), A for both these two 
patterns has been plotted as the average of the pairing 
amplitudes on interacting and non-interacting sites. Due 
to the proximity effect, even in the absence of interaction 
on a lattice site, there exists a finite value of pairing am- 
plitude through the tunneling effect from its neighboring 
interacting sites. While Ay for both the checkerboard 
and stripe patterns consistently increases as a function 
of U, A (u=o) increases up to C/(r=i) ~ 2.0 where accord- 
ing to Fig. ^a), the homogeneous and inhomogeneous 
lattices have the same pairing amplitude. A(u=o) then 
falls off for larger U values. Hence, the region of growth 
of A(u=o) with U coincides with pairing amplitude of the 
inhomogeneous system being large. 

The superconducting transition temperature of a lat- 
tice with a bimodal Ui can also be larger than a uni- 
form interaction distribution. According to Fig. Q^a), 
a,t U — 1.5, A(T — 0) is enhanced due to inhomogene- 
ity but not at U — 3.0. Fig. ^ illustrates the collapse 
of A as the temperature is increased for the uniform, 
checkerboard and stripe patterns. In Fig. EJa) corre- 
sponding to U = 1.5, inhomogeneity has a similar ef- 
fect on the gap through all values of finite T up to Tc 
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which is significantly larger for the inhomogeneous pat- 
terns. In Fig. ^h) corresponding to t7 = 3.0, at T = 0, 
A for the uniform pattern is significantly larger than 
its inhomogeneous counterparts. However, because the 
gap of the inhomogeneous system is nonzero to higher 
temperatures than the uniform system, for values of 
0.5 ^ T ^ 1.25, there exists an enhancement region 
for A for inhomogeneous patterns. For the uniform pat- 
tern, we find the Tc in very good agreement with the 
BCS value, k^^T^ « (A(0)?7)/1.76, as expected for our 
mean field treatment. Note that the effect of inhomo- 
geneity on Tc would appear to be even more dramatic 
if we compared the inhomogeneous curves against their 
BCS counterparts with the same A(r = 0). A similar 
increase in Tc upon introducing a checkerboard pattern 
is found in the MCMF calculations as well, arising from 
the loss of long-range phase coherence. This is partic- 
ularly significant because the MCMF incorporates the 
subtle nature of the superconducting transition in 2D 
discussed earlier. We have also independently confirmed 
that our conclusions and arguments equally apply for a 
model with nearest-neighbor attraction, leading to a d- 
wave SC close to half-filling, which reflects the cuprates' 
phenomenology more truthfuUyj^ 

IV. SUMMARY AND DISCUSSION 

In summary, we have shown that for the attractive 
Hubbard model on a square lattice, there is a significant 
range of doping and interaction strength over which the 
superconducting order parameter is larger for Ui = 0, 2U 
than for Ui = U uniformly. It is worth emphasizing that 
in most situations, inhomogeneities reduce values of order 
parameters and critical temperatures, even when com- 
parisons are made, as they are in this paper, to homoge- 
neous systems with the same average value of all param- 
eters. This is true, for example, of classical site diluted 
Ising models, where the ferromagnetic J is increased to 
compensate for absent sites, and quantum models like 
the boson Hubbard model where random chemical po- 
tentials monotonically decrease and ultimately destroy 
superfluidityi^SiSi An exception is the increase of TnccI 
by randomness reported in DMFT studies of the repul- 
sive mo del 1^ 

The increase in the SC gap has been verified for the 
checkerboard, stripe and random patterns and thus is 



insensitive to the pattern of disorder. The growth is 
due to the proximity effect, i.e., the tunneling effect of 
the Cooper pairs from the interacting sites leading to fi- 
nite order parameter values even on non-interacting sites. 
This conclusion is supported by the effect occurring at 
weak coupling, where the coherence length is large, rather 
than in the strong coupling regime of preformed pairs. 
Agreement between the BdG results and the MCMF 
calculations justifies the application and conclusions of 
the BdG approach within the small C7 regime. Finally, 
nonuniform interaction strength also leads to the strong 
signal in the phase transition temperature Tc- Counter- 
intuitively, this increase in Tc continues even for values of 
U for which the the order parameter is larger with uni- 
form Ui than a bimodal choice. However, in this weak 
coupling parameter regime, Tc is a supralinearly increas- 
ing function of [/, so it may be that in the inhomogeneous 
system, the sites with larger U produce a nonlinear en- 
hancement relative to Tc of the homogeneous system and, 
through the proximity effect, drag the Ui — sites along 
with them. 

While the attractive Hubbard Hamiltonian obviously 
does not incorporate many of the features of high Tc su- 
perconductors (notably the symmetry of the pairing), the 
model has been shown to provide useful insight into some 
of their phenomenology, for example the spin-gap. It is 
therefore tempting to speculate that our results concern- 
ing inhomogeneity may have similar connections. Specifi- 
cally, recent ARPES data^ suggests that the underdoped 
phase of LSCO (La2-a:Sra:Cu04) consists of SC clusters, 
embedded in the AF host. In such a system, inhomoge- 
neous gaps appear naturally and our results here indicate 
that the superconducting transition is in fact determined 
by the largest gap values rather than the much smaller 
gaps found at phase boundaries, as one would naively 
think. This renders the SC phase more stable than it 
would otherwise be, and also simplifies the description of 
these systems. 
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